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ABSTRACT 

A rational function is presented for the acoustic spectra generated 
by deflection of engine exhaust jets for under-the-wing and over-the- 
wing versions of externally blown flaps. The functional representation 
is intended to provide a means for compact storage of data and for data 
analysis. The expressions are based on Fourier transform functions for 
the Strouhal normalized pressure spectral density, and on a correction 
for reflection effects based on the N-independent-source model of 
P. Thomas extended by use >&f a reflected ray transfer function. Curve fit 
comparisons are presented for blown flap data taken from turbofan engine 
tests and from large scale cold-flow model tests. Application of the 
rational function to scrubbing noise theory is also indicated. 


INTRODUCTION 

Future short takeoff and landing (STOL) aircraft will undoubtedly 
be required to meet stringent noise level regulations. Thus, noise is 
one of the primary factors in the selection of a STOL propulsion-lift 
system. One candidate source of propulsive -lift is the externally blown 
flap (EBF) system. Propulsive lift is generated by downward deflection 
of engine exhaust jets in either the tmder-the-wing or the over-the-wing 
version of EBF. Each EBF STOL concept introduces noise sources not 
present in conventional takeoff and landing (CTOL) airplanes which cause 
an increase in noise. These EBF STOL propulsive lift concepts are being 
investigated by scale model tests (refs. 1 to 8), and full scale engine 
tests (refs . 9 and 10) . Concurrently, theoretical and experimental 
studies are being conducted on various noise source mechanisms (refs. 11 
to 14), and empirical prediction methods are being devised (ref. 15). 

The objective of this paper is to determine a rational function 
representation of flap noise spectra that is useful (1) as an empirical 
equation to fit EBF STOL spectral data, and (2) as a function that can 
be related to certain noise-source mechanisms. The rational function 
representation is intended to serve as a link between experimental EBF 
spectral data and theoretical explanations, and to provide a mathematical 
expression for EBF noise spectra that can become part of an empirical 
prediction method. 
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Many noise source mechanisms may contribute to the overall noise of 
an EBF system. These include: (1) leading edge noise caused by incident 

turbulence; (2) scrubbing noise generated in the turbulent boundary layer 
produced by the jet mixing region and convected along the wing or flap; 

(3) separated flow noise; (4) trailing-edge noise caused by vortex 
shedding of eddies from the deflected flap; (5) jet mixing noise from the 
distorted deflected exhaust jet; and (6) engine internal and turbomachin- 
ery noise. Each of these sources may dominate at some angle depending on 
the specific jet exhaust velocity and EBF STOL configuration. 

The noise source mechanism receiving most emphasis is the one to 
which scrubbing noise is attributed (refs* 16 through 25). In spite of 
many advances in understanding this and other noise sources, procedures 
for calculating noise from EBF STOL airplanes starting with basic aero- 
acoustic mechanisms are not available currently. Hence, present design 
studies depend on model and full scale tests and on empirical relations 
determined from these tests. 

For design studies, one of the most important aspects of the meas- 
ured noise is the third-octave sound pressure level spectrum designated 
in this report by SPL. It is used to compute the perceived noise level 
(PNL) and the effective perceived noise level (EPNL) , which are the 
present criteria used to judge the effect of an aircraft on the community 
noise level. 

This report describes an empirical rational function that can be 
used to represent the SPL data. The function provides a representation 
for broadband SPL spectra. With the representation described, sound 
spectra can be stored in a computer in a compact form and comparison of 
large amounts of data from different experiments is facilitated. The 
representation can also be used to extrapolate data to untested condi- 
tions, which is useful in making estimates of PNL, EPNL, and footprint 
calculations. In functional form, new relations or physical insights 
into source mechanisms may more easily become apparent. In fact, the 
empirical functional representation is related herein to a theoretical 
description of scrubbing noise, and thus provides a reason for the ob- 
served similarity of EBF data. In achieving this functional representa- 
tion, corrections for acoustic reflections in ground test data are taken 
into account. 

The paper will present the development of equations representing the 
spectra. A basic equation is given in terms of the spectral density of a 
real random process. This equation will be related to scrubbing noise 
theory. Next, a simplified equation that represents the spectra will be 
given. Equations for corrections for reflection and atmospheric effects 
will also be included. The complete calculation method will then be dis- 
cussed, and finally, results from EBF data taken from turbofan engine 
tests and cold flow tests of under- the-wing and over- the -wing test models 
will be given. 
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ANALYSIS 


Approach 

The approach taken to represent flap noise spectra is based upon 
finding a function of Strouhal number, such that the third— octave 

sound pressure level (SPL) can be fitted by 

SPL(i) = 10 (1) 

where 

OASPL = 10 logj^Qjsj^ j^qSPL(JO/10J ^2) 

D is an effective diameter (refs* 7, 8, and 15) 
e 

V is an effective velocity (refs* 7, 8, and 15) 
e 

and 

i refers to the third octave frequency band 


(All symbols are defined in the symbol list.) Specifically, 


^ ^ j27T 



j2TTSt 


(3) 


The requirements on the function <J>(i2) are: 

(1) It should be tailored for applicability to large Strouhal num- 
bers* Curve fits that are accurate at large Strouhal numbers are neces- 
sary, since externally blown flap engines are likely to have large ef- 
fective diameters and low exhaust velocities. 


(2) It should be consistent with theoretical considerations. 

(3) It should depend on as few mathematical coefficients and opera- 
tions as possible since this will reduce computer calculation time. 

(4) It should be applicable to a wide range of system configurations 
(e*g*, single or coannular nozzles, mixer nozzles, over-the-wing, and 
under-the-wing, etc*) . 

From equation (1) , 

10 log^Q{((>(n)} = SPL(£) - OASPL - 10 logj^QjAf^(Jl) ~ 




(4a) 
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The quantity represented by equation (4a) is the Strouhal normalized mean- 
square pressure spectral density level per unit Strouhal number. In this 
paper expression (4a) will be referred to as the PSD“*" function, PSD"*" 
which is a function of Strouhal number is calculated from SPL data free 
of atmospheric absorption and reflection effects (lossless data) , It has 
been shown that spectral data taken over a range of exhaust velocities 
for numerous engine-over-the-wing and engine -under- the -wing externally 
blown flap configurations correlate when plotted in terms of the FSB"*" 
function (e,g,, refs, 3 and 4), Thus, 

PSD'^ = 10 log^Q[(f,(J2)] (4b) 

where cj)(^^) is seen to be the Strouhal normalized mean square pressure 
spectral density per unit Strouhal number. 

Since the acoustic noise generation is basically a random process, 
a function form for c{)(Q) will be selected such that it represents the 
power spectral density of a real random process. The expression will 
also be capable of fitting a wide range of PSD"^ data. This basic func- 
tional form will be related to a theoretical description of scrubbing 
noise. Next, the basic expression will be simplified by removing certain 
factors. The simplified expression will also fit the PSD"^ curve, but 
having fewer terms it can be evaluated faster in a computer calculation. 

Finally, the equations necessary to relate the calculated lossless 
PSD'^ curve to SPL data affected by ground reflections and atmospheric 
attenuation will be stated. The basis for the correction is discussed, 
and some examples of the ground reflection effect are presented. 


Spectral Density Function 

Basic function . - The mathematical criteria that c|)(n) must satisfy 
to represent a power spectral density is that it must not be negative and 
must be integrable. However, to represent the power spectral density of 
a real random process, (f)(Q) must also be an even function of Strouhal 
number. A basic function having these properties can be forinulated from 
the reciprocal of the Fourier transform of a ’’damped cosine oscillation.” 
Explicitly, the function is 


o(^2,c) 


- c) (Q - c*) 

Q + 


(5) 


where c represents a complex number and the asterisk (*) indicates the 
conjugate relation. 


The function a(Q,c) defines a curve of a for a given value of c. 
Since the desired expression for <()(^^) involves the logarithmic addition 
of a number of individual o(Q,c) curves (to fit observed PSD"^ variation). 
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the general expression for a real random process, <|)Real^^^ » 
pressed as the ratio of the products of these terms, 




■ ^ 

a TT a(£J,aj^) 

° i*l 


I ^ I ® *^lr ^ 

k=l 


( 6 ) 


a is a normalization constant; the complex numbers aj^^ and bj^ have 
negative real values and positive or zero imaginary values; and is 

less than . The individual cuirves comprising the total function 
4iReal(^) are referred to as Stroujial response functions. The unknown in 
equation (6) is the number of response functions and their value (a^^ and 
bj^ terms) required to properly fit a given data variation. 

Application to scrubbing noise theory . - Equation (6) can be used to 
yield inputs to a theoretical description of scrubbing noise. Starting 
from Curie's extension of Lighthill’s theory to surfaces (ref. 16), it is 
possible to show that dipole broadband noise is radiated from a rigid 
surface adjacent to a turbulent boundary layer. The noise produced is 
referred to as scrubbing noise. The mean square pressure per unit 
Strouhal number from references 24 and 25 would correspond to 


(}>{J2) = 



ip(£5)A^(SJ)dS(x) 


^—00 

^ / St^ i|»(n)A^(£5)dS(x)dSt 


(7) 



-^- T ~ 

in the notation of this paper, where S(x) represents the surface im- 
mersed in a flow, is the spectrum of the pseudo-sound or fluctuating 

surface pressure, and A^.(n) is the correlation area over which the nor- 
malized cross spectral: density of the normal stress is unity. 


Some attempts have been made to define equations for and for 

A(.(fi) based on experimental data for surfaces other than flaps (refs. 24 
arid 25) . The necessary measurements do not exist specifically for flap 
noise sources. However, empirical representations of the ifi(fl) and 
Aj,(fi) functions for flap noise have been made based on experimental data 
of similar phenomena. In reference 14, Fink approximated the correlation 
area (in the notation of this paper) by 


A^,(i^) 



2ir 


(£2 - a ^) 


|2 


(8a) 
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and the static pressure fluctuation spectrum was approximated by 




2ir(U - a^) 

(Q - bi) (fi - b2) 


(8b) 


where 


bf = irSt^^ (-1 - j) 
= TrSt^/2 (-1 + j) 


Expressions for these two functions can also be postulated from 
equation (6) as follows: 






i=l 


a(f2,aj^) 


K. 


I I 


m=J 


(9a) 


and 




K(j (2 

7T 


m=l 


Nb 

I I 0(^5 ,bjj^) 


k=l 


(9b) 


where ^ and has a negative real part and a zero or pos- 

itive imaginary part. 

The expressions of equation (9) are more general and provide greater 
flexibility than the expressions of equation (8) , The greater number of 
coefficients in equation (9) compared to equation (8) increase the poten- 
tial ability of the functions to fit experimentally based and 

functions and should also be more useful in fitting PSD data. 

Simplified function , - Equation (6) is the power spectral density of 
a real random process. A modification of equation (6) that yields a 
simpler expression which is equally suitable for representing PSD"*" data 
will be discussed next. However, while the new equation can be used to 
satisfactorily fit a spectral density it does not represent a power spec- 
tral density of a real random process. 
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The modification consists of dropping the term - c )/(f^ + / cc^) 
from equation (5) . The equation for the calculated lossless PSD*^ then 
can be written as 


PSD 


+ 

cal 


= 10 logj^Q 



2 


af) 


N. 




- b^) 


( 10 ) 


where all symbols are the same as in equation (6) . 


Reflection Correction 

Acoustic reflections cause problems > since in some experiments re- 
flections distort the true sound spectrum and lead to incorrect spectral 
shapes with frequent over estimation of the overall sound pressure level 
(OASPL) . The reflection correction developed herein is based on a mathe- 
matical model. Reflection model parameters are adjusted simultaneously 
with the functional representation parameters of equation (10) by the 
computer pr'ogram. The calculated SPL depends on both sets of parameters. 
In situations where the acoustic reflection model parameters are unknown, 
the functional representation plays a vital part in correcting reflection 
effects. A description of the reflection model and some examples are 
provided herein. 

The measured sound pressure level, ^thout atmospheric a.ttenua-“ - 
tioip obtained by including reflection effects with the lossless value of 
SPL is given by 


SPL(f^) = SPL'^(f^) + 10 |T(f^)|^ (11) 

where |T(fp)|^ is the reflectance and SPL"*"(f^.) (eq. (1)) is the sound 
pressure level without atmospheric attenuation and ground reflection. 

In arriving at equation (11) the model assumes the noise source is 
distributed over a region and represents this by rays from N independent 
sources. This assumption was used by P. Thomas (ref, 26) to study acous- 
tic interference of the noise produced by a single jet due to reflections 
from a plane. An alternative treatment for a distributed source using 
the wave equation is given in reference 27. The ground reflection effect 
on a single point noise source has been treated in references 28 to 31. 

The working hypothesis will be that the noise can be treated as if 
it were from a vertical arrangement of N independent sources of equal 



8 


strength at the trailing edge of the last flap. This model permits one 
to consider a microphone at a particular angle to be receiving reflec- 
tions from the wing, flap and ground surfaces, resulting in an angle- 
dependent interference effect. 

The sound reaching the microphone is assumed to consist of the sound 
that travels the direct path, going a distance sj_, and the sound that 
travels a reflected path going a distance S 2 (fig. 1) . Since the re- 
flected path is longer than the direct path, the reflected sound reaches 
the microphone after a time delay x where x = (s2 - s ^) /cq relative 
to the direct sound that left at the same time. Thus if the direct wave 
reaching the microphone was emitted at a time t, the reflected wave 
reaching the microphone at the same time is due to one emitted at a time 

t - X . 


Many point sources are assumed to be causing the sound and each 
source is assumed to produce reflections from wing, flap, and ground sur- 
faces. This causes many reflection paths. The amplitude of the reflected 
sound is changed by these multiple reflections before the signal reaches 
the microphone. The modification to the reflected pressure signal will 
be assumed to be related to the direct ray by a convolution of the rer- 
flecting surfaces' impulse response with the direct pressure signal: 


P'(t) 4 



g(x)p(t - x)dx 


( 12 ) 


The Laplace transform of the reflecting surface impulse response 
g(x) is the reflected ray transfer function G(s) . This complex re- 
flected ray transfer function can also be written in terms of its phase 
6 and amplitude [q[. Thus, 


G(s) = !q(s)|e^'^^®^ k 



(13) 


The source, flap, and microphone geometry with respect to the ground 
plane is shown schematically in figure 1. Each of the N-independent 
sources will be assumed to produce a direct ray that travels a distance 
®1* All the sources will be assumed to be at a fixed distance r from 
the microphones. For the n^li source at a height hg(n) 

r? 9^3■/2 

Si(n) = {r^ + [hg(n) - h^]^} 

The reflected ray path distance for the n^h source is 


S2(n) = (r2 + [hg(n) + 
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The observed relative pressure at a distance r will be given by 
the sum of the pressure signals traveling directly from the N sources 
and the reflected pressure signals which are delayed by a time delay 
relative to the direct signals- and reduced by a divisor, z, proportional 
to the ratio of the travel distance of the reflected pressure signal to 
the direct pressure signal. The variation of the heights of the source, 
relative to the radial distance is assumed small; thus the mean value of 
the ratio of the reflected signal travel distance to the direct signal 
travel distance will be used. 


The resulting mean-square pressure spectrum density at the micro- 
phone is found by a two-step procedure. First the auto-correlation of 
the signal is found. Then the Laplace transform of the auto-correlation 
is taken. This two-step procedure is used to obtain the frequency re- 
sponse since the noise will be assumed to be produced by a stationary 
Independent random process. 


The spectrum without attenuation is assumed to be the sum of the 
spectra of the N sources. These sources are also assumed to be of 
equal strength. If the source spectrum, , |(5(s) | and 6(s) are assumed to 
be constant over any third-octave band, the thirdi-octave spectrum (ne- 
glecting atmospheric attenuation) measured at the microphone is given by 
equation (11) where the reflectance, [T(f(,)P, is given by 


|T(f^) 


= 1 + 


lQ(fJ 


2 lQ(fc>l ^ 

" ^ L 

n-1 


N 


cos[27rf X 
. c n 




sin(7i Af T ) 
c n 

(14) 


For use in the determination of |x(f^)|^ by a computer, the re- 
flected ray transfer ftmction G(s) will be assumed to be independent of 
frequency. Hence in equation (14) , | Q(fcH will be replaced 

by the constants |q| and 6. 


Typical reflectances for a low and high level of |q| are shown in 
figures 2 and 3. Figure 2 applies to the test configuration at the cold 
flow model test site (ref, 4) and figure 3 applies to the test configura- 
tion at the EBF test site (refs, 9 and 10), Table I contains the values 
of the constants used in the reflectance calculations. The effect of low 
|q| values is to minimize the magnitiide of the cancellations and rein- 
forcements. 


CALCULATION METHOD 

The complete equation used to represent the SPL at a microphone loca- 
tion for a given configuration is based on equations (1), (4), (10), (11), 
and (14), 
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OASPLp/10 D 
10 ■ Af 

c V 

e 



X 1 + + |- ' COs(27Tf T -- 6) 


8111(71 Af T ) 
c n 

7T Af T 
c n 


where A<j«(f^,r) is the atmospheric attenuation. The atmospheric attenua- 
tion was computed by the method of reference 32. 

A computer program is used to adjust the parameters OASPLp, a^, a^; 
i “ 1, • . b]j^; k = 1, . • IqI and 6. A diagram of the 

method is shown in figure 4. The parameters are adjusted to minimize the 
total ”cost,” C^, consisting of the sum of the squares of the differences 
between the measured and calculated SPL values, C^, and a penalty cost 
Cp (refs, 33 and 34) due to selection of parameters that violate the 
parameter constraints previously mentioned. (See eq. (2) for 0ASPL+ con- 
straint and eq. (6) for and b|^ constraint.) The total cost func- 

tion is nonlinear. The problem of minimizing the nonlinear cost function 
was solved by search techniques. These techniques consist of systematic 
procedures for varying the cost function parameters until a minimum value 
of the cost function is found. The basic methods available are discussed 
in references 33 and 34- To provide flexibility in the choice of 
and and due to the presence of parameter constraints, a search tech- 

nique which does not require evaluation of derivatives was chosen. The 
search technique was that of Powell (ref. 35). The computer program used 
was adapted from reference 36. The optimization procedure was previously 
applied successfully to a similar identification problem in reference 37. 

The input information for the method consists first of the constants 
required for the reflection model. These consist of the following: the 

number of sources, N, which was selected to be 5; the mean value of the 
ratio of the distance the reflected ray travels to the distance the di- 
rect ray travels, z; and the time delay, associated with each of the 
five sources. Second, the constants and % for the PSoJ^^i model 

are needed. For the data considered was selected to be 3 and 
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was selected to be 1. Last, initial values of the parameters to be varied 
(OASPlJ, a^, |q|, 6, a^, bg) are necessary. 

RESULTS 

The calculation method based on the rational PSD’*’ function with re- 
flection corrections was applied to several sets of measured data to il- 
lustrate the comparisons obtained for the PSD"^ function and for both meas- 
ured and corrected SPL spectra. Also included are examples of OASPL var- 
iations with exhaust velocity and PSD“^ functions over a range of micro- 
phone angles . 


Basic PSD^ Function 

Several sets of lossless blown flap data are available for determin- 
ing fits by the PSD"*" function based on the formulation for a real random 
process given by equation (6) . Data correlations over a range of exhaust 
velocities in terms of lossless j}pectral density versus Strouhal 
number are available in reference 3 for a cold flow test of an under— tha- 
wing externally blown flap configuration. Comparable Strouhal correla- 
tions for a cold flow over-the-wing configuration with attached flap flow 
(powered lift) and with unattached flow (conventional lift) are presented 
in reference 4. Calculated values of PSL"** for these data were obtained 
from iterative calculations of equations (6) and (4b) in conjunction with 
the measured values of PSD"*". A total of nine parameters were used in the 
calculation. The parameter values corresponding to the best fit of the 
data are given in table II. 

Figure 5 shows a comparison of the data and the curve fit for PSD^ 
for the under -the -wing model of a flap setting of 30*^ ^ 60^ for the 85^ 
microphone angle position. Figure 6 shows a comparable comparison at a 
microphone angle of 100^ above an over-the-wing model with attached flow 
at a flap setting of 10^ - 20° ♦ The comparison for the unattached flow 
version of the over-the-wing model with flap settings of 10° - 20° and 
30° - 60*^ are shown in figure 7« In all cases, the calculate curve 
fits are in very good agreement with the data over the entire Strouhal 
number range. Thus, equation (6) should be applicable to data from many 
different types of blown flap configurations. 


Simplified PSD Function 

The simplified PSD"*" function of equation (10) was also used to fit 
the data shown in figures 5 to 7. The nine parameters obtained are shown 
in table III in terms of the complex constants needed for the fit. In 
all cases, the curve fit to the data was very nearly as good as for the 
basic function of equation (6) . An example of the comparison is given in 
figure 8 which shows the calculated PSD“*" values for the basic (eq. (6)) 
and simplified (eq. (10)) functions as obtained for the fit of the data 
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of figure 6. On the basis of the close correspondence of the two re- 
sults, it was concluded that the simplified version of the PSD"^ function 
given by equation (10) is completely adequate for routine calculation 
purposes . 


Nutther of Parameters 

The curve fits used in fitting the data shown in figures 5 to 7 were 
based on using = 1 and % = 3, t^ich results in nine parameters in 
equations (6) and (10) . Figure 9 is a plot of the individual Strouhal 
response functions F(aj^), F(bj^), F(b 2 ), F(b 3 ) corresponding to a^^, bj^, 
b 2 , and b 3 that make up the simplified PSD+ functions of figure 8. 

The b 3 curve is seen not to contribute to the final curve form. In 
fact, if the magnitude of b^ (i = 1, 2, or 3) is larger than 60 and the 
maximum Strouhal nunher of interest is less than 50, the value of F(b^) 
i s nea rly zero and a(J^,b£) is nearly constant and has a value of 
>^ib^. Therefore, b^^ can be eliminated and the number of parameters re- 
quired can be reduced from nine to seven. Tables II and III show that 
for each of the four curves, the magnitude of one of the b^^ is larger 
than 60. Consequently a seven -parameter formulation of equation (6) or 
(10) should be adequate for many calculations. 


Application to Scrubbing Noise Theory 

Two sets of functions were presented earlier for the correlation 
area A<.(fi) and the static pressure fluctuation spectrum i)>(J^) for equa- 
tion (7) describing flap surface scrubbing noise. One functional form of 
<J)(£2) based on using equation (8a) for Aj,(fJ) and equation (8b) for i|/(fi) 
with Stg = 0.3 was developed by Fink (ref. 14). The particular form of 
the denominator parameters selected by Fink results in a (j)(f2) function 
that depends on only a single constant, St^; however, in the notation of 
this paper, this functional form consists of three parameters. In arriv- 
ing at this total the normalization constant, aQ is counted as one param- 
eter of (j)(n) . Also the functional form of <|)(fi) will be said to have 
two parameters in the denominator. This is done since a complex constant 
and its complex conjugate are counted only once and equation (8) indicates 
the two parameters (bj^ and b2) in the denominator of <)>(i^) are conjugate 
complex constants. 

The other functional form based on equations 9(a) for A^(fl) and 
9(b) for ip(fi) is the one developed in this paper. This functional form, 
as has been previously discussed, was found to give a good fit to the 
data with = 1 and N^ = 2, to yield a total of seven parameters. 

Thus, the major characteristic of equation (9) compared to equation (8) 
is the addition of two parameters to the numerator and two parameters to 
the denominator with freedom to select the six best parameters for the 
numerator and denominator. 

In figure 10 the ability of two different functional forms of iji(n) 
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to fit the data shown in figure 5 is compared. Each PSD"^ curve shown in 
figure 10 is calculated from equation (4) which relates the PSD curve 
to the ^RealW function and from equation (7) which relates the 
'^’Real^^^ function to the correlation area, and the static pressure 

fluctuation spectrum • The three-parameter functional form of <|>(^^) 

derived by Fink with St^ set equal to 0.3 was found to be in general 
agreement with the data of figure 5 but overestimates the PSD"^ at large 
Strouhal numbers as mentioned by Fink (ref. 14). This can be seen in 
figure 10 where the seven-parameter curve that provides an exact fit to 
the data of figure 5 is compared with three-parameter curve fit. Thus, 
the greater number of coefficients in the functional form of 4^Real^^) 
developed in this paper and the use of the computer to find the best 
parameters provides an excellent curve fit to PSD"^ data. 


Application to Measured Data 

Cold flow model . - Measured and calculated values are compared in 
figure 11 for a number of calculated quantities for the case of data from 
an over-the-wing model (ref. 4) with two wing flaps and attached flap 
flow (powered lift configuration). The test configuration is for cold 
flow with a circular nozzle and a flap setting of 10 ^ 20 ^# 

Figure 11(a) shows the measured third-octave sound pressure level 
spectrum. The data were taken at an exhaust velocity of 253 m/sec 
(830 ft/sec) and a microphone angle of 100^ from the inlet. The cancel- 
lations denoted by C]^, 02,03 and the reinforcements denoted by R 2 ,R 3 »R 4 
are due to reflection effects. Figure 11(b) shows the measured sound 
pressure level, SPLj^^, the calculated SPL including effects of atmospheric 
absorption and reflections, and the calculated lossless sound 

pressure level, reflectance used to account for the reflec- 

tion effect is that shown in figure 2(b). Comparison between measured 
(SPL^) and calculated (SPL^^l) spectra indicates that SPLq^i gives a 
good fit to the data. 

Figure 12 illustrates the fit between calculated values of lossless 
spectral density, and the PSD"*" spectrum obtained from the meas- 

ured SPL by using the calculated reflectance shown in figure 2(b) and a 
calculated atmospheric attenuation. This calculated PSD based on the 
measured SPL data is denoted as FSdJ. The close proximity of the least 
square curve fit, PSdJo^]^, to the PSD^ calculated points is apparent. 

Full-scale engine configuration . - Another example of the applica- 
tion of the method is given for a set of data from the tests of refer- 
ence 9 for a full-scale EBF system with a three-flap wing and a turbo- 
fan engine with a coannular nozzle. The test condition is for a takeoff 
flap setting of 0^-20^-40^, an effective exhaust velocity of 242.3 m/sec 
(795 ft/sec) (maximum power), and a microphone angle of 80° from the in- 
let. Figure 13(a) shows the measured third -octave SPL spectrum for the 
test condition. Reflection effects are indicated by the first cancella- 
tion denoted by and by the second reinforcement denoted by R 2 . The 
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turbofan engine in this program was highly noise suppressed (ref. 10) 
with much of the fan dxscrete tone noxse removed. Howeverj because some 
residual fan tone noise may be present at frequencies above 4000 hertz 
for the maximum power operating test condition only data from 50 to 
3150 hertz were analyzed herein. 


The calculated spectrum for lossless sound pressure level, SPL"*" 
is shown in figure 13(b). Also shown are calculated values of SPL 'in- 
cluding reflection and atmospheric attenuation effects, SPL^al, and values 
of the measured sound pressure spectrum, SPL^. The reflectance used to 
account for the reflection effects is similar to that shown in fig- 
ure 3(a) . Comparison of SPLcal with the measured SPL data, SPLjj, 
indicates a good representation of the data was achieved. 

Figure 13(c) illustrates the fit between PSdJ and PSD^gi. A good 
least-squares curve fit was again achieved. 


Variations with angle and velocity . - Calculations were also made 
with the iteration method developed using data from the full-scale engine 
EBF configuration previously discussed taken at microphone angles from 
20 to 140° over a range of exhaust velocities. Figure 14 shows values 
of OASPLp determined by the curve fitting procedure. The lossless 
values of OASPL were calculated from the value of lossless SPL. Slopes 
of 0ASPL+ following a sixth or an eighth power of the effective ve- 
locity are shown in the figure for comparison purposes. 

According to figure 14, for microphone angles between 20° and 90° 
the OASPLp varies with the sixth power of the effective exhaust velocity 
reflecting the dominance of surface scrubbing (or dipole) noise. For 
angles between 100° and 140°, the data in figure 14 indicate a trend 
toward an eighth power dependence on effective velocity. At 120° the 
eighth power relation fits the data quite well. This is interpreted to 
be an Indication that the jet (or quadruple mixing) noise source is more 
xmportant at these greater angles. 


Figure 15 shows 




for angles from 40° to 120° at the maximum 


o ''' 

power setitng wxth = 242 m/sec (795 ft/sec) . The curves for micro- 
phone angles of 40 to 80° (sixth power of velocity region) are similar 
xn appearance and peak at a Strouhal number between 0.25 and 0.27. The 
^®°cal curves for microphone angles of 100° and 120° (eighth power of 

^®Sion) have their peak values at a Stroiahal number between 0.30 
and 0.35. All the PSD+ curves fall within a band of less than 4 dB 
width. 


CONCLUDING REMARKS 

A calculation procedure based on a simple rational function for 
spectral density has been shown to provide a representation for the 
acoustic spectra generated by deflection of engine exhaust jets for under- 
e-wing and over-the-wing versions of externally blown flap configura- 
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tlons. The representation is useful in correcting reflection effects, 
and provides a compact method of storing sound spectral data. In addi' 
tion, the simple function is related to a more basic pressure spectral 
density function that can also be used to fit PSD data and can be 
related to scrubbing noise theory. 

The simple rational function combined with the procedure for cor- 
recting reflection effects yields mean-square pressure level per unit 
Strouhal number curves that are devoid of interference effects and can 
be easily compared with one another. 


SYMBOLS 






a . 

X 


C 

P 


c 


c 

o 


D 


m 

F(c) 


f 


G(s) 

g(^^) 


correlation area over which the normalized cross spectral 
density of the normal stress is unity 

atmospheric attenuation, dB 

zero of i = 1, 2, . • • > 

a 

a complex number, normalization constant 
pole of k = 1, 2, , * ^ 

cost function, sum of square of error between measured and 
calculated sound spectra 

penalty function 

total cost function 

cancellation frequencies 

a complex number 

velocity of sound, m/sec 

effective diameter, m 

a complex constant, m=l, 2, . * *>^d 

Strouhal response function 

frequency. Hertz 

nominal mean frequency, Hertz 

reflected ray transfer function 

reflecting surfaces^ impulse response function 
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h 

o 

hg(n) 


S, 

% 

N 




microphone height, m 
th 

height of n source, m 

number of a(fi,d^) terms in A^(Q) and t);(n) 
third-octave band number 

number of numerator o(fi,a. ) terms of di , (n) 

1 Real 

number of independent noise sources 

number of denominator a(£2,b,) terms of ij> - (^2) 

source index 


OASPL 

oaspl'*' 

P 

PSD 

PSO**" 




p 

P’ 



r 


S(x) 


overall sound pressure level relative to p dB 

lossless OASPL parametric value, dB 

Strouhal normalized mean-square pressure level per unit 
Strouhal number, dB 

PSD without atmospheric absorption and ground reflection 
effects, dB 

PSD"*^ calculated from formula, dB 

PSE"*” calculated from measured data, dB 

PSD calculated from a real random process spectrum 
density, dB 

2 

pressure, N/m 

2 

reflected pressure signal, N/m 

reference pressure 2x10“^ N/m^ (2x10"^ microbar) 
magnitude of reflected ray transfer function 
reinforcement frequency 
distance from source to microphone 
surface immersed in a flow 


SPL 

spl”^ 


third-octave sound pressure level reference to p dB 

third-octave sound pressure level referenced to p that 
would be measured without atmospheric attenuation and in 
the absence of reflecting surface, dB 
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SPL 


cal 


SPL. 


M 


St 


St 

o 

s 




V 

e 

X 


z 


6(s) 

0 

a(fl, c) 


T 

T 

n 




SPL calculated from an equation, dB 
SPL measured, dB 
Strouhal number, 

constant, Strouhal number (0.3, ref. 14) 
j27Tf, jo) 

distance from source to microphone along direct ray path, m 
distance from source to microphone along reflected ray path, m 
reflectance 
time , sec 

effective velocity, m/sec 
surface position vector 

mean ratio of path length reflected ray traveled to path 
length direct ray traveled for all n sources 

real constant, -27rSt^ 

bandwidth at center frequency f (0.2316 f ), hertz 

c o 

reflected ray transfer function phase factor, rad (deg) 
microphone angle relative to inlet 

function of Strouhal number such that the reciprocal of the 
Fourier transform of a "damped oscillation" is proportional 
to I a(Q,c) I ^ 

time delay, sec 

time delay between arrival of wave leaving source at time t 
and traveling direct path to receiver and one leaving at 
time t and traveling reflected ray path to receiver, sec 

Strouhal normalized mean-square pressure per unit Strouhal 
number 




Strouhal normalized mean-square pressure per unit Strouhal 
number for a real random process 




spectrum of the pseudo-sound or fluctuating surface pressure 
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JJ relative angular velocity, j27rf(D^/Vg) = j2irSt 
(jj angular velocity, 2-irf, rad/sec 


n 

product of terms with index i ranging from m to n 

i=m 

n 

^ ^ summation, sum of terms of index i ranging from m to n 
l«m 

Superscript: 

* denotes conjugate 
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TABLE I, - CONSTANTS USED IN REFLECTANCE CALCULATIONS 


(a) Constants used for reflectance shown in figure 2. 
N = 5; r = 15.24 m (50 ft); h,,. = 3.88 m 
(12.75 ft) 


n 


1 

2 

3 

4 

5 

hs(n) 

m 

4.130 

4.0081 

3.886 

3.764 

3.642 


ft 

13.55 

13.15 

12.75 

12.35 

11.95 

"n 

sec 

0.0058 

0.00567 

0.0055 

0.0053 

0.0052 


(b) Constants sued for reflectance shown in figur 
N = 5; r = 30.48 m (100 ft); h^ = 12.743 m 

e 3. 

(9 ft) 

n 


1 

2 

3 

4 

5 

hs(n) 

m 

3.475 

3.736 

2.743 

2.377 

2.012 


ft 

11.4 

10.2 

9.0 

7.8 

6.6 

^n 

sec 

0.00182 

0.00164 0 

.00145 

0.00125 

0.00106 



TABLE II. - PARAMETERS TO CALCULATE SELECTED PSD^gal USING EQUATIONS (6) AND (4b) 

Dimensionless Complex number, C (real value, imaginary value), C = a^, a,, b, , b«, or fa- 

parameter 

Data of fig. 5 Data of fig. 6 Data of fig. 7(a) Data of fig. 7(b) 

(151,610, 151.610) (133.3510, 133.3510) (372.2495, 372.2495) (293.9065, 293.9065) 

a^ (-0.1965826, 0.1824718) (-0.1304851, 0.1214311) (-0.7096163, 0.6316123) (-0.4733036, 0.2092519) 

(76.22998, 26.7030) (-0.1540204, 0.4026451) (-0.3045551, 1.036240) (-0.4142318, 0.7848299) 

(-0.9577132, 0.7647091) (-1.564607, 1.282075) (-3.892744, 3.885970) (-4.280862, 3.136002) 

b^ (-1.546595, 2.192013) (-76.433486, 35.98859) (-156.9277, 15.68585) (-87.00853, 66.38647) 

CO 

(F- 

TABLE III. - PARAMETERS TO CALCULATE SELECTED PSD^^ CURVES USING EQUATION (10) 

Dimensionless Complex number, C (real value, imaginary value), C = a-, a-., b, , b«, or b-, 

parameter 

Data of fig. 5 Data of fig. 6 Data of fig, 7(a) Data of fig. 7(b) 

(192.6902, 192.6902) (105.7856, 105,7856) (290.5190, 290.5198) (268.0784, 268.0784) 

a^ (-0.2528992, 0.1505393) (-0.1947860,0.06866252) (-0.1126150,0.08250151) (-0.1289076, 0.0709157) 

b^ (-1.194895, 0.5249010) (-0.1479430, 0.4224659) (-0.4334988, 0.3890225) (-0.4899988, 0.424448) 

b^ (-1.745644, 1.715307) (-1.963396, 1.474141) (-3.177522, 0.5692431) (-3.651598, 0.3990379) 

b^ (-115.7714, 14,24272) (-73.48344, 45.69005) (-124.6561, 34.93369) (-98.45943, 52.17599) 


PRECEDING NOT FILMED 




Figure 1. - Schematic of source and microphone geometry 
relative to ground. 



(bl loi'1.0, 6-5.6. 

Figure 2. - Calculated reflectance ieq. 114)) for test site used to 
obtain cold flow blown flap data in reference 4 evaluated for a 
low and high level of |Q I . 
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Figure 3. - Calculated reflectance (eq. (14)) for test site 
used to obtain turbofan engine external ly-blown -flap 
data (ref. 9) evaluated for a low and high level of |Q|. 







Figure 4 - Diagram of computer method. 
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Figure 5. - Comparison of curve fit using equation (6) with 
Strouhal correlation of cold flow under the-wing externally- 
flap data taken from reference 3. Flap setting, 30^ - 
microphone angle, 85^, 
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Figure 6. - Comparison of curve fit using equation i6) with 
Strouhal correlation of data for engine-over -the-wing 
externally -blown -flap configuration with powered lift taken 
from reference 4. Flap setting 10° - 20° microphone angle, 
100 °, 

/ ' • 


EXHAUST VELOCITY, 



(b) FU\P POSITION. 30° *60° 


Figure?, - Comparison of curve fit using equation (6) 
with Strouhal correlation of data for the engine-over- 
the wing externally-blown -flap configuration with con- 
ventional lift taken from reference 4. Microphone 
angle, 120° 



Figure 8. - Comparison of use of basic and simplified 
spectral density functions to fit data of figure 6. 






Figure 9. - Strouhal response functions used to produce 
simplified PSO'*' function in figure 8. 




(b) CALCULATED AND MEASURED SOUND PRESSURE LEVEL 
SPECTRUM. 

Figure IL - Over-the-wing cold-flow powered lift model with flap 
position 10° - 20° {ref. 41. Exhaust velxity, 253 m/sec 
1830 ft/sec); distance, 15.2 m (50 ft). Microphone height, 3.9 m 
<12.75 ft); source height, 3, 9 m {12.75 ft); microphone angle, 
100°. 



Figure 10. - Plot of PSD'*' functions based on data shown in 
figure 5. Comparison is made between three parameter 
curve fit based on equation (8) and seven parameter curve 
fit based on equation (9), 



Figure 12. - Comparison of measured and calculated spectral 
density using over-the-wing cold-flow powered lift model 
data of figure 11. Vg - 253 m/sec (830 ft(sec). 
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(b) MEASURED AND CALCULATED SOUND PRESSURE 
LfVEL SPECTRUM. 

Figure 13. - Turbofan engine externa I ly-blown-f lap 
configuration with effective exhaust velocity, 

242 , 3 m/sec (795 ft/sec); distance, 30, 48 m (100 ft); 
microphone and source height, 2,74 m (9 ft); 
microphone angle, 80^. 



STROUHAL NUMBER, St, HD^/Vg) 
(c)PSD+g, COMPARED WITH PSDJ^. 



EFFECTIVE EXHAUST VELOCITY, Vg, M/SEC 
(g) 0 = 100°. (h) 0 = 120° 0)0 = 140°. 


Figure 14, - Variation of lossless overall sound pressure level 
with effective velocity at microphone angles, 0, from 20° to 
140° for turbofan engine EBF configurations. Flap setting. 0°- 
20°'40° Microphone radius, 30.48 m (100 ft). 


Figure 13. - Concluded. 
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Figure 15, - Variation of calculated lossless spectral 
density PSDggi at effective exhaust velocity of 
242.3 m/sec u95 ft/sec) with microphone position 
over a range of angles from 40^ to 120®. Turbofan 
engine EBF configuration, flap setting, 0®~20®-40®. 


